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February 1, 2008 

Abstract. The Nagel-Schreckenberg model is a simple cellular automaton for a realistic description of 
single-lane traffic on highways. For the case Vmax ~ 1 the properties of the stationary state can be obtained 
exactly. For the more relevant case Vmax > 1, however, one has to rely on Monte Carlo simulations or 
approximative methods. Here we study several analytical approximations and compare with the results 
of computer simulations. The role of the braking parameter p is emphasized. It is shown how the local 
structure of the stationary state depends on the value of p. This is done by combining the results of 
computer simulations with those of the approximative methods. 

PACS. 02.50.Ey Stochastic processes - 05.60.-|-w Transport processes; theory - 89.40.+k Transportation 



1 Introduction 

Cellular automata (CA) do not only serve as simple model 
systems for the investigation of problems in statistical me- 
chanics, but they also have numerous applications to 'real' 
problems pi. Therefore it is not surprising that in recent 
years CA have become quite popular for the simulation of 
traffic flow (see e.g. [|[^]). 

CA are - by design - ideal for large-scale computer 
simulations. On the other hand, analytical approaches for 
the description of CA are notoriously difficult. This is 
mainly due to the discreteness and the use of a parallel 
updating scheme (which introduces 'non-locality' into the 
dynamics). In addition, these models are defined through 
dynamical rules (e.g. transition probabilities) and usually 
one does not have a 'Hamiltonian' description. Therefore 
standard methods are not applicable. Furthermore, one 
has to deal with systems which do not satisfy the detailed 
balance condition. 

However, there is a need for exact solutions or, at least, 
for good approximations. These results as well as other 
exact statements may help to greatly reduce the need for 
computer resources. The interpretation of simulation data 
is often difficult because of the 'numerical noise' and finite- 
size effects. Even in the cases where an exact solution is 
not possible, a combination of analytical and numerical 
methods might provide better insights. This is especially 
true for non-equilibrium systems where only a few exact 
or general results exist which could serve as a guideline. 

In recent years, several analytical approximation meth- 
ods for the Nagel-Schreckenberg model |4] for single-lane 
highway traffic have been proposed J5l@, 7,0. All of these 
approximations yield the exact result for the stationary 
state in the case Vmax = 1- However, these investiga- 
tions focused on the so-called fundamental diagram, i.e. 
the flow-density relation. Here we will reinvestigate these 



approximations and calculate further quantities of inter- 
est in order to determine the accuracy. We will also dis- 
cuss several limits and emphasize the effect of the brak- 
ing probability p. We focus on the effect of p on the mi- 
croscopic structure of the stationary state. Such micro- 
scopic characterizations have recently been used success- 
fully for the asymmetric exclusion process |^,|0|. Apart 
from Vmax — 1 here only the case of Vmax = 2 is inves- 
tigated since one does not expect a qualitively different 
behaviour for Vmax > 2 ||6[- 

The paper is organized as follows: First we will briefly 
recall the definition of the Nagel-Schreckenberg model in 
Sec. y and the different analytical approximations in Sec. 
0. In Sec. several physical quantities, e.g. the funda- 
mental diagram, headway and jam-size distributions and 
correlation lengths, are calculated using the analytical re- 
sults. In Sec. 1^ the predictions of the approximations are 
compared with each other and with results from computer 
simulations. The local structure of the stationary state is 
investigated as a function of the randomization p. In the 
final Sec. ra a summary of the results together with our 
conclusion are given. 



2 The Nagel-Schreckenberg model 

The Nagel-Schreckenberg (NaSch) model H] is a proba- 
bilistic cellular automaton. Space and time are discrete 
and hence also the velocities. The road is divided into 
cells. The length of a cell is determined by the front- 
bumper to front-bumper distance of cars in the densest 
jam and is usually taken to be 7.5 m. Each cell can either 
be empty or occupied by just one car. The state of car 
j (j = 1, . . . , N) is characterised by an internal parame- 
ter Vj (vj = 0, 1, . . . , Vmax), the instautaueous velocity of 
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the vehicle. In order to obtain the state of the system at 
time i + 1 from the state at time t, one has to apply the 
following four rules to all cars at the same time (parallel 
dynamics) B: 



+ l,w„ 



.) 



Rl Acceleration: Vj{t + 1/3) = mm{v j{t) 

R2 Braking: Vj{t + 2/3) = mm{dj{t),v-i{t + 1/3)) 

R3 Randomization: Vj{t + 1)— max(wj(i + 2/3) — 1,0) 

with probability p 
R4 Driving: car j moves Vj{t + 1) cells. 

Here dj (t) denotes the number of empty cells in front of car 
j, i.e. the so-called headway. For Vmax = 5 a calibration of 
the model shows that each timestep t ^ t+1 corresponds 
to approximately 1 sec in real time B. For simplicity we 



special variants of the NaSch model, e.g. the VDR model 
[ [13| in the slow-to- start limit. 

In the following the analytical approaches are discussed 
briefly for the cases Vmax = 1 and Vmax — 2. All meth- 
ods other than MF theory are exact for Vmax = 1- For 
Vmax = 2, on the other hand, they are only approxima- 
tions. 

Applications of the analytical methods to the calcula- 
tion of physical quantities, e.g. the fundamental diagram, 
jam-size distributions and correlation lengths, are given in 
Sec. |. 



3.1 Mean-Field Theory 



1 







1 




2 







7,5 meter 



Fig. 1. Configuration in the Nagel-Schreckenberg model. The 
number in the upper right corner gives the velocity of the car. 



will consider only periodic boundary conditions so that 
the number of cars is conserved. The maximum velocity 
Vmax can be interpreted as a speed limit and is therefore 
taken to be identical for all cars. Fig. [^ shows a typical 
configuration. Throughout the paper we will assume that 
the cars move from left to right. 

The four steps have simple interpretations. Step Rl 
means that every driver wants to drive as fast as possible 
or allowed. Step R2 avoids crashes between the vehicles. 
The randomization step R3 takes into account several ef- 
fects, e.g. road conditions (e.g. slope, weather) or psycho- 
logical effects (e.g. velocity fluctuations in free traffic). An 
important consequence of this step is the introduction of 
overreactions at braking which are crucial for the occu- 
rance of spontaneous jam formation. Finally, step R4 is 
the actual motion of the vehicles. 

The NaSch model is a minimal model in the sense that 
all four steps R1-R4 are necessary to reproduce the basic 
properties of real traffic. For more complex situations (e.g. 
2-lane traffic [Q or city traffic |12 ) additional rules have 
to be formulated. 



3 Analytical Methods 

In this section several analytical approaches which have 
been used for the description of the NaSch model are 
reviewed. The simplest, a mean-field (MF) theory, com- 
pletely neglects correlations. Since MF theory turned out 
to be inadequate even for the fundamental diagram, im- 
proved methods have been developed which allow to take 
into account short-range correlations exactly. All methods 
described here are microscopic theories since macroscopic 
theories are not able to describe the NaSch model prop- 
erly. However, they are extremely useful and accurate for 



The simplest analytical approach to the NaSch model is 
a (microscopic) mean-field (MF) theory [|| . Here one con- 
siders the density Cy (j, t) of cars with velocity v at site 
J and time t. In the MF approach, correlations between 
sites are completely neglected. 

For Vmax = 1 the MF equations for the stationary state 
{t^oo) read §: 



Co = {c + pd)c , 
ci — pcd 



(1) 
(2) 



with c = Co -\- ci, d = 1 — c and p — 1 — p. The flow is 
simply given by fMF{c) = ci. 

For random-sequential dynamic^ the MF approach is 



known to be exact for u„ 



1 



For parallel dynamics. 



however, MF theory underestimates the flow considerably 
(see Sec. |1|). 

For Vmax — 2 the rate equations for the densities are 
given by |6) 



Co = (c -I- pd)co + 
ci = d [pcQ + {pc - 
C2 =pd^{ci -I-C2) 



(1 +pd)c{ci -f C2) 
Vpd){ci +C2)] , 



(3) 
(4) 
(5) 



The solution is given by 



Co = 



Cl 



C2 



i-pd^ ' 

p{l — pd'^)dc 
1-pd^ 
p'^d^c 



(6) 



1-pd^ ' 



and the flux can be calculated using fmric) = ci + 2c2. 
Again the flow is underestimated considerably (see Sec. 
4.1). This is true for arbitrary Vmax- The general form of 
the MF equations can be found in S. 



^ In random-sequential dynamics in each timestep a cell 
which is updated is picked at random. 
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3.2 Cluster Approximation 



with the conditional probabiHties 

P{ri 



following we will use occupation variables tj where tj =0, 
if cell j is empty, and tj — v, if cell j is occupied by a car 
with velocity v. The change of the ordering has to be taken 
into account in the calculation of observables. The flow is 
given by /(c) =pP(l,0). 

For an n— cluster of n consecutive cells with state vari- 
ables T^" (t) = {tj (i), . . . , Tj+n-i (i)) at time t it is straight- 
forward to derive an exact evolution equation for the clus- 
ter probability P{Tj{t), . . . ,Tj+n-i{t))- One has to take 
into account that cars can enter the cluster from one of the 
Vmax cells to the left of the cluster and can leave the clus- 
cells to the right. Therefore the state 



ter to one of the Vmax 



(n) 



at 



. ^ (t+l) of a cluster depends not only on its state t^ 
time t, but also on the neighbouring cells iTj-v„^ax (J')i ■ ■ ■ i 
Tj-i{t)) and (Tj+„(i), . . . , r,-+„+«„„^_i(t)). 

In the stationary state the master equation for clus- 
ter r^"^ = [tj^--- ,'^i+n-i) of n cells has the following 
structure: 



(7) 



r(n + 2u„a^) 



where x("+2^™'") = (tj-«„,, , ■ 
sition probabilities V7(t("+^^" 



■ >'^j+n+u„„,-i)- Thetran- 
= ) __> x^"^) have to be de- 
termined from the rules R1-R4. Note that the master 
equation for n— clusters involves (n -I- 2i)maa;)— clusters. If 
the stationary state is translation-invariant the probabili- 
ties P{tj . . . ,Tj^n-i) are independent of j. 

In order to obtain a closed set of equations one has to 
express the (n -I- 2umaa;)— clusters through the n— cluster 
probabilities. At this point some approximation has to 
be made. Usual one uses a factorization into products of 



-clusters. We illustrate this for v„ 



Fig. I): 



2 and n = 3 (see 



The cluster approximation |g|,|6| is a systematic improve- 
ment of MF theory which takes into account short-ranged 
correlations between the cells. In the n-cluster approxi- 
mation a cluster of n neighbouring cells is treated exactly. 
The cluster is then coupled to the rest of the system in a 
self-consistent way. Related approximations have already 
been used (under different names) for other models M, 
|l|,0,|l7). 

In order to simplify the description, we sometimes choose 
a slightly different update-ordering R2-R3-R4-R1 instead 
of R1-R2-R3-R4, i.e. we look at the system after the ac- 
celeration step. Then there are no cars with v — and 
effectively we have to deal with one equation less. In the 



P{ti\ t2,.. . ,T„ 

P{ri, 



J ■ • ■ I In/ 



Li ■ ■ ■ J 



P{Tl,...,Tn) 



(9) 
(10) 



If we denote the probability to find the system in a 
configuration (ri, . . . , r^) by P(ti, . . . , tl) the 1-cluster 
approximation means a simple factorization 



(11) 



P{n,...,rL) = l[Pir,). 

This is nothing but the mean-field theory of section ^.l[ 
For the 2-cluster approximation one has a factorization of 
the form 

P(ti,... ,tl) (X P{ti,T2)P{t2,T3) ■ ■ ■ P{tl-utl)P{tl,ti) . 

(12) 

The 3-cluster approximation is depicted graphically in 
Fig. |c. 



b) 












Fig. 2. Graphical representation of the n— cluster approxima- 
tion for a) n = 1 (i.e. mean-field theory), b) n = 2, and c) 
n — 3. Shown is the central n— cluster starting at site j and 
all clusters which have to be taken into account in the master 
equation (Q) for Vmax = 2. 

In general, the master equation in n-cluster approx- 
imation leads to {vmax + 1)" nonlinear equations. This 
number can be reduced by using the so-called Kolmogorov 
consistency conditions [|5[ 

^ P(ti,... ,T„_i,r) = P(ri,... ,r„_i) 



r=0 



E^(-' 



Tl, 



,T„„i) . (13) 



P{t^'^) = P(r,-2 h-i,r,- ) . P(r,-i| T„T,+i ) • 

• -P(Ti>Tj+i.Tj+2) • (8) 

• ■P(tj+i,Tj+2|tj+3) • P(rj+2,rj+3|Tj+4) 



Nevertheless, a solution is only feasible for relatively small 
cluster-sizes [pi. The quality of the approximation im- 
proves with increasing n and for n —t oo the rt-clustcr 
result becomes exact. 
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p = 0.5 




Fig. 3. Fundamental diagram for Vmax = 1: Comparison of 
computer simulations (•) with the exact solution (full line) 
and the mean-field result (broken line). 
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Fig. 4. Comparison of simulation results with the n-cluster 
approximation {n — 1, . . . ,5 from bottom to top) for the fun- 
damental diagram for Vmax = 2 and p = 1/2. 



However, for Vmax = 1 already the 2-cluster approx- 
imation is exact MM- The 2-cluster probablities for the 
stationary state are given exphcitly by 



P(0,0) = l-c-P(l,0), 
P(l,l) = c-P(l,0), 



(14) 



P(1,0) 



^(0'^) = ^ 



1- Vl - 4pc(l - c) 



where againp = 1—p. The flow is given by /(c) = pP{l, 0). 

Higher order cluster approximations with n > 2 yield 
the same result (H). This indicates that (O) is exact. 
Indeed this has been proven in H by a combinatorial ar- 
gument. 

For Vmax = 2 already the 2-cluster approximation yields 
a nonlinear system of equations which has to be solved nu- 
merically. The fundamental diagrams obtained from the 
n-cluster approximation (n = 1,... ,5) are compared in 
Fig. with results of Monte Carlo simulations. One can 
see a rapid convergence and already for n ~ 4 the dif- 
ference between the simulation and the cluster result is 
extremely small. 



3.3 Car-Oriented Mean-Field Theory 

The car-oriented mean-field (COMF) theory H is another 
possibility to take into account correlations in an analyt- 
ical description^]. 

The rules R1-R4 can be rewritten in terms of dj and 
Vj only. Therefore the state of the system can be charac- 
terized instead of the occupation numbers tj equivalently 
by the headways and the velociticfQ The central quanti- 
tiy in COMF is the probability Pn(v) to find exactly n 
empty cells (i.e. a gap of size n) in front of a car with 
velocity v. In this way certain longer-ranged correlations 
are already taken into account. The essence of COMF is 
now to neglect correlations between the headways. 

Using again the update ordering R2-R3-R4-R1 (see 
Sec. 3^), one obtains for Vmax = 1 the following system 
of equations resulting from the master equation for the 
stationary state (with P„ = P„(u = 1)): 



Po=g[Po+pPi] , 

Pi = gPo + [Pg + pg] Pi + PgP2 , 

Pn = pgPn-1 + [pg + pg] Pn + PgPn+1 



(15) 



{n> 2) 



where g = pJ2n>i ^" ~ Pi^ ~ ^o] is the probability that 
a car moves in the next timestep. g — 1 — g then is the 
probability that a car does not move. 

As an example for the derivation of these equations 
we consider the equations for n > 2. Since the velocity 
difference of two cars is at most 1, a gap of n cells at 
time t + I must have evolved form a gap of length n — 
1, n, or n -|- 1 in the previous timestep. A headway of 
n — 1 cells evolves into a headway of n cells only if the 
first car moves (with probability g) and the second car 
brakes in the randomization step (probability p), i.e. the 
total probability for this process is pgPn-i- Similarly, the 
headway remains constant only if either both cars move 
(probability p(7) or both cars do not move (probability p^). 
Finally, the headway is reduced by one, if only the second 
car moves (probability p^). 

The probabilities P„ have to satisy the normalization 
conditions 



- oo 



(16) 
(17) 



n=0 



(O) is a consequence of the conservation of the number of 
cars since each car with headway n occupies n + 1 cells. 

Although the infinite system of non-linear equations 
( [l5| ) looks more difhcult than those of the cluster approxi- 
mation, a solution possible using generating functions |f7|. 



A similar method for reaction-diffusion models is discussed 



m]l 

3 



For a complete equivalence one has to track the position 
Xi of one car, e.g. car 1. From xi, {dj} and {vj} one can then 
determine {tj}. 
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For Vr, 



1 one finds 



Pn = 



P„ 



1 
2pc 

p {Po+pil-Po) 



2pc - 1 + Vl - 4pc(l - c) 
P(l - Po) 



(n>l), (18) 



and for the flow f{c,p) — eg again the exact solution is 
reproduced. 

For Vmax — 2 one has two coupled systems of the type 
(p^), since one has to distinguish P„(i; — 1) and Pn{v = 2) 
and the probabilities ga that a car moves a — 0,1,2 cells 
in the next timestep. This system has a similar structure 
as (^5|) and is given explicitely in 0. It can also be solved, 
but does not give the exact solution. 

Note that the COMF approach assumes an infinite sys- 
tem size. For a finite system of length L the largest head- 
way that can appear is M := L — N . For Vmax = 1 it 
is possible to satisfy ( [l6| ) and ( jl7| ) by just changing the 
equations for Pm-i and Pm'- 



Pm-1 = PgPM-2 

Pm ^pgPM-i 



[P9+P9]Pm~i + 
[l-pg+ pg] Pm 



P9]Pi 



M 



(19) 



The equations for Pn with n < M — 1 are the same as 
in (^5|). In order to derive (|lt^ ) we have made the Ansatz 
Pm-1 = aiPM-2 + a2PM-i+ofiPM and Pm = f3iPM-i + 
/32Pm and determined the cofficients aj and j3j from (Hq) 



nations for all 



and (|1_ 

For Vmax = 2 it appears that the 
have to be modified in order to satisfy (|16| ) and ( |17[) . This 
corresponds nicely to the qualitative difference of the cases 
Vmax = 1 and Vmax > 2 and the ap pear ance of long-ranged 
correlations in the latter (see Sec. 3.4). 



3.4 Garden of Eden States 

An important effect of the parallel dynamics is the exis- 
tence of configurations which can not be reached dynam- 
ically [g|. These states are called Garden of Eden (GoE) 
states or paradisical states since one never gets back once 
one has left |1^. An example for a GoE state is given in 
Fig. ||. Note that the velocity is equal to the number of 
cells that the car moved in the previous timestep. For the 
configuration shown in Fig. ^this implies that the two cars 
must have occupied the same cell before the last timestep. 
Since this is forbidden in the NaSch model, the configura- 
tion shown can never be generated by the dynamics. 




30 




Vmax =2; p = 


= 0.5 












- 


25 
20 


: 


Simu .-, 


*v\ 




- 


15 


: 


i /-'MF 


"■"■■■■^ 




'-_ 


10 


- 


/' 






-__ 


05 


/ 




III 


1 


% 



Fig. 6. Comparison of tlie fundamental diagrams obtain from 
MC simulations with the results of mean-field theory (MF) and 
mean-field theory without GoE states (paradisical mean-field, 
pMF). 



GoE states. One can therefore hope that by eliminating 
all GoE states and applying mean-field theory in the re- 
duced configuration space (paradisical mean-field, pMF) 
one will find a considerable improvement of the MF re- 
sults. 

For Vmax = 1 all states containing the local configura- 
tions (0, 1) or (1, 1) are GoE states, i.e. the cell behind a 
car with velocity 1 must be emptyn. This only affects eq. 

and the equations for pMF theory read: 



m 



Co = 7V(co -I- pd)c , 
ci — Afpcd , 



(20) 
(21) 



where the normalization A/" ensures cq + ci = c and is given 
exphcitly by A/" = l/(co -|- d). 

Solving (|2|) for Ci by using cq = c — ci , one obtains ci 
as a function of the density c: 



ci = 1 (l - v/1 - 4p(l 



(22) 



(with p = 1 — p). Since the flow is given by /(c) = 
we recover the exact solution for the case Vmax — 1 



= Ci 

(see 



Section 3.2). 

It is interesting to note that for random-sequential dy- 
namics MF theory is exact whereas for parallel dynamics 
pMF is exact. Therefore the origin of the correlations is 
the parallel update procedure. The existence of the GoE 
states is responsible for the differences between parallel 
and random-sequential dynamics. This is probably not 
only true for the NaSch model, but is a rather general 
property of CA models. 

For Vmax — 2 one has to take into account more GoE 
states [||. pMF is no longer exact, but it still leads to a 
considerable improvement of the MF results (see Fig. H). 



Fig. 5. A Garden of Eden state for the model with Vmax > 2. 



The simple mean-field theory presented in the previ- 
ous section does not take into account the existence of 



* Here we do not use the changed update order, so '0' denotes 
a cell occupied by a car with velocity 0. 
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4 "Physical" Quantities 

In this section we will use the analytical methods of Sec. || 
to determine physical quantities of the NaSch model, e.g. 
the fundamental diagram, cluster-size distributions and 
correlation functions. These results will be used in Sec. ^ 
to gain a better understanding of the microscopic struc- 
ture of the stationary state. 



4.1 Fundamental Diagram 

For Vmax = 1 the fundamental diagram in MF approxi- 
mation is given by 



fMF{c) =pc(l 



(23) 



As already mentioned in this case 2-cluster approximation, 
COMF and pMF are exact and the flow is given 



4.2 Headway Distribution 

Fig. shows the result for the distribution of headways, 
Pn, for Vmax — 1- Sincc in this case COMF is exact, we 
can use the result dlq). The headway distribution has just 
one maximum, located at n = 1 for small densities and at 
n = for large c. 

These two regimes also exist for higher velocities. At 
high densities the headway distribution is maximal a,t n = 
0, whereas for low densities the maximum is found at some 
value Umax > 0, where Umax depends on the density. In 
addition, an intermediate density regime exists, where dis- 
tribution exhibits two local maxima, one at n = 0, cor- 
responding to jammed cars, and one at nioc > 0, corre- 
sponding to free flowing cars pO| . 



Vmax =5; p=o.5 



1 



/(c,p) = -(l-v/l-4p(l 



c)c 



(24) 



(with p = 1 — p). The fundamental diagram is symmet- ^ 
ric with respect to c ~ 1/2, reflecting the particle-hole 
symmetry of the NaSch model with v,nax = 1- 

It is apparent from Fig. 0that MF considerably under- 
estimates the flow for intermediate densities. This shows 
that correlations are important in this regime. These cor- 
relations lead to an increase of the flow. One finds a particle- 
hole attraction (particle-particle repulsion), i.e. the prob- 
ability to find an empty site in front of a car is enhanced 
compared to a completely random configuration. In terms 
of cluster probabilities this particle-hole attraction can be 
expressed more quantitatively as P(1,0) > P(1)P(0) = 
c(l-c). 
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As mentioned in Sec. 3.4, pMF yields the exact solu- 
tion for Vmax — 1- This shows that there are no 'true' 
correlations apart from those due to the existence of GoE 
states. The 2-cluster approximation and COMF then be- 
come exact since both methods are able to identify all 
GoE states. 

For Vmax > 1 one expects stronger correlations than 
for Vmax = 1 §• Indeed for Vmax = 2 the difference be- 
tween the MF prediction 



fMF{c) = Ci +2C2 



p{l+p(P)dc 
l~p(P 



(25) 



and the Monte Carlo data is even larger (see Fig. Q). How- 
ever, pMF is not exact and therefore 'true' correlations ex- 
ist. This is in agreement with the fact that MF is also not 
exact for the NaSch model with Vmax > 1 and random- 
sequential dynamics. 

The correlations can be taken into account systemati- 
cally by the n— cluster approach which converges rapidly 
with increasing n. Already for n = A one obtains a very 
good agreement with the simulation results (Fig. 0). The 
quality of COMF depends strongly on the value of p. In 
general its prediction is between the 2- and 3-cluster ap- 
proximation. 



Fig. 7. Distribution of headways for different densities for 
Vmax = 1 and p = 0.5 (left) and Vmax = 5 and p = 0.5 (right). 



4.3 Correlation Length 

As an application of the cluster approximation we will 
compute density-density correlations for Vmax = 1 in the 
following. Using occupation numbers Uj =0,1 the density- 
density correlation function is defined by 

(mn,.) = 2^ ninrP{ni, ... ,71^) (26) 

{rij} 

where the prime indicates that the sum runs over all states 
with fixed particle number TV = m -t- • • • + n^ . 

In order to evaluate the sum in (pq), it is convenient 
to use a grand-canonical description. One introduces a 
fugazity z which controls the average number of particles 
{rij) and sums over all configurations in (P6|). Using the 
2-cluster approximation (O) - which is exact for Vmax — 1 
- the correlation function is given by 



{nirir) 



— ^niUrZ^ P{ni,n2)P{n2,n3) ■ 



{".} 



■ P{nL-i,nL)P{nL,ni) (27) 
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with N = J2i 



and the normahzation 



4.4 Jam-Size Distribution 



Zgc = X! ^^ n ^K'"j+l)- 



(28) 



Introducing the transfer matrix 



f p(o,o) V^P(o,i) 

lyiF(l,0) zF(l,l) 



this can be written succinctly as 
Zg, = Tr P^ , 



{ninrj 



^TrfqP'-iQP^-'^+i) 



(29) 

(30) 
(31) 



with the matrix Q{ni,n2) = niP(ni,n2). 

The correlation length ^ can be obtained from the 
asymptotic behaviour (r — ^ oo) of the correlation func- 
tion 



{niUr) — c oc e '"'^ 



(32) 



where c = (rij) is the (average) density of cars. ^ is de- 
termined by the ratio of the eigenvalues X± of P (with 
|A+|>|A_|): 



r' = ha 



A- 



A, 



The explicit expression for £_ can be obtained from 
1 



^-=2 



A±y/A^+4pz[P{l,0)]^ 
with A = P(0,0) + zF(l,0) . 



(33) 



(34) 
(35) 



P{a,b) are the cluster probabilities given in (14). The 
fugazity z can be related to the density c via the equa- 
tion c = {rij). Using an expression analogous to (|3l| ) one 
finds 



cA+ = 



B{[P(l,0)]2 + i3P(l,l)} 



[P(l,0)]2 
with P = A+ - P(0, 0) . 



L52 



(36) 

(37) 



For fixed p, f is maximal for c = 1/2 which corre- 
sponds to z = 1. ^(c — 1/2) diverges only for p ^ 0. In 
that case one finds that ^(c = 1/2) oc p~^/^. Monte Carlo 
simulations show that the same behaviour still occurs (at 
density c = l/{vmax + 1)) for v^ax > 1 ||2^. Therefore, 
the correlation function already gives an indication that 
the system is not critical for p > 0. The simulations show 
that there is no qualitative difference between the cases 
Vrnn^ — 1 and Vmax > 1 as far as the phase transition is 



concerned |21 



The above results demonstrate that, although the 2- 
cluster approximation is exact, not all correlation func- 
tions are of finite range. 



For a better understanding of the differences between the 
cluster approximation and COMF we calculate the distri- 
bution of jam-sizes using these methods. Let C„ be the 
probability to find a (compact) jam of length n, i.e. n 
consecutive occupied cells. 

In COMF, C„ is proportional to (1 - Po)Po • • • Po(l - 
Po) = (1 — Po)^Pq^^^. The number of jams is proportional 
to 1 — Pq so that one finds 



L'n 



(1 - Po)PS- 



(38) 



This distribution is purely exponential and P„ > P„+i for 
all n. Therefore COMF is not able to describe clustering 
or phase separation, i.e. situations where jams with more 
than one car dominate. Clustering implies that there are 
correlations between gaps which are completely neglected 
in COMF. 

In 2-cluster approximation, C„ is proportional to P(0| 1) 
P(l|l)--- P(l|l)P(l|l)P(ilp) where we have used the 
conditional probabilities (E^). and the different update- 
ordering R2-R3-R4-R1. Note that the case n — 1 has to be 
treated separately, Ci ex X]II=r -f(0|^)^fe|0)- The num- 
ber of jams is obviously given by Mj — Yl^v'=i P{v\^) ^^d 
one obtains 

C(^^ = ;^ P(0|1)P(1|1)"-^P(1|1)P(1|0) {n > 2), 



ci'^^^f^pimpim 



Mj 



(39) 



.(2) 



Cn decays exponentially for n > 2, especially one has 
Cn+i l£ Cn- For the ?7i-cluster approximation one can de- 
rive similar expressions. Since here clusters of size m are 
treated exactly it is now possible to find the maximum of 
the jam-size distribution somewhere between n = 1 and 
n = m. For n > m the jam-size probability decays expo- 
nentially due to the mean-field-like self-consistent coupling 
of the cluster to the rest of the system. 



4.5 Further Applications 

We briefly mention other results which have been obtained 
analytically using the cluster approximation. In |2y] the 
distributions of gaps and the distance between jams have 
been calculated for Vmax = 1 using the 2-cluster approx- 
imation. For the gap distribution one recovers the (ex- 
result which 



act) result which has been derived in Sect. 3.3 (see eq. 
(Iq)). For the calculation of the distribution of gaps be- 
tween jams one defines every vehicle with velocity to be 
jammed. The distance between two jams is then given by 
the distance between a car with velocity and the next 
car with velocity 0. In |p3 the probability V{t) of a time 
headway t has been investigated. This quantity is defined 
in analogy to measurements on real traffic where a de- 
tector registers the time interval between the passing of 
consecutive cars. For a discussion of the behaviour of the 
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time-headway and the jam-gap distribution in the NaSch 
model we refer to |2y,^. These quantities are also very 
useful for studying variants of the NaSch model which ex- 
hibit phase separation P3[. 



5 Microscopic Structure of the Stationary 
State 

In the following we demonstrate that the microscopic struc- 
ture of the stationary state changes qualitatively with the 
braking parameter p. For illustration we first discuss the 
deterministic limits p — and p — \. After that the be- 
haviour for randomizations < p < 1 is discussed. How- 
ever, we concentrate on the limits p <C 1 and 1 — p <C 1. 
Finally, in Sec. 5.4 we compare COMF and cluster ap- 
proximation. 



Here we focus on the ability to reproduce 
the correct microscopic structure of the stationary state. 



5.1 p = 

In the absence of random decelerations the velocity of the 
vehicles is determined solely by its headway. For densities 
c < l/(wma2; + l) all headways can be larger than Vmax and 
therefore the cars move with velocity Vmax in the station- 
ary state. This is no longer possible for c > l/{vmax + !)• 
Here the velocity of the cars is determined by the average 
headway d available, d = [L - N)/N = (1 - c)/c. The 
fundamental diagram is therefore given by 



f{c) = 



VmaxC for C < l/{Vmax + 1), 
1-C iOT C>l/{Vmax + l)■ 



m 



Note that in general the stationary state is not unique, 
but determined completely by the initial condition. 

For p = there is no tendency towards clustering |gj| . 
Overreactions are not possible and there is no spontaneous 
formation of jams. In this sense the case p = is unrealis- 
tic. The dynamics is completely determined by "geomet- 
rical" effects, since the behaviour of the cars only depends 
on the available average headway. There is no 'attractive' 
interaction between the cars and therefore no mechanism 
for clustering. 



5.2 p = l 

In the case p — 1 the dynamics is again deterministic. 
However, it is very different from that of a model with 
maximum velocity v^ax — 1- The point that we want to 
stress here is the existence of metastable states. For Vmax = 
2 and density c = 1/3 the state ..1..1..1.. (here '.' denotes 
an empty cell and '1' a cell occupied by a car with velocity 
1) is stationary with flow /(c = 1/3) = 1/3. On the other 
hand, the state ..0..0..0.. is also stationary with vanish- 
ing fiow, since a standing car will never start to move for 
p=l. 

For densities c > 1/3 all stationary states have vanish- 
ing flow. Here at least one car has only one empty cell in 



front. Therefore after step R2 this car has velocity 1 and 
will then decelerate to velocity in step R3. After that it 
will never move again. 

For densities c < 1/3 stationary states with non-zero 
flow exist. These are not stable under local perturbations, 
i.e. stopping just one car leads to a complete breakdown 
of the flow. In this sense these states are metastable. 

Starting from random initial conditions, the metastable 
states have a vanishing weight in the thermodynamic limit, 
since already one standing car leads to a zero-flow state. 

Note that the metastable states at p = 1 do not exist 
for Vmax = 1. In that sense the difference between Vmax = 
1 and Vmax > 1 becomes most pronounced at p = 1 . 

The deceleration introduces a kind of 'attractive' inter- 
action between the cars which can lead to the formation 
of jams. However, these jams are typically not compact, 
but of the form .0.0.0.0.. A car approaching a standing 
car adapts its velocity in step R2 such that it reaches the 
cell just behind the standing vehicle. In step R3 it then 
decelerates further and so there will be a gap of 1 between 
the cars. 



5.3 <p< 1 

In the limit p — > 0, the fundamental diagrams obtained 
from COMF and the 3-cluster approximation become ex- 
act, in contrast to the 2-cluster approximation (see Fig. ph. 
Even for values of p « 0.1 there is an excellent agreement 
between the fundamental diagrams obtained analytically 
and the Monte Carlo simulations. Since 'realistic' values 
of p are in the region p ~ 0.1 — 0.2, the approximations 
are indeed applicable in the relevant parameter regime. 





Fig. 8. Fundamental diagram for Vmax = 2 and p — 0.1: Com- 
parison of simulation data (•) with the 2-cluster (• • ■ ) and 3- 
cluster results ( — ) (left) and with COMF ( — ) (right). 



Fig. H shows a comparison of the computer simulations 
with the COMF results for the probabilities ga that a car 
moves a sites. Although there is an excellent agreement 
for the fundamental diagram /(c) = c{gi + 2g2) at p = 0.1 
(see Fig. @), there are considerable deviations at interme- 
diate densities for the individual values of gi and 32- This 
is somewhat surprising and indicates that COMF is not 
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exact in the limit p — > 0, i.e. that there are stiU correla- 
tions b etwe en the headways. This will be discussed further 
in Sec. 



5.4 





Fig. 9. Comparison of the COMF prediction (• ■ ■ ) with com- 
puter simulations ( — ) for the probabilities ga that a car moves 
a sites, go can be obtained from go — 1 — gi — g2- The braking 
parameter is p = 0.1 and v^ax = 2. 



On the other hand, the 3-cluster probabilities obtained 
from our simulations show a very good agreement with the 
results of the 3-cluster approximation (Fig. 10[). This sug- 
gests that the 3-cluster approximation is asymptotically 
exact for p ^ 0. 

The limit p ^ 1 is difficult to investigate numeri- 
cally. Preliminary results show that COMF is not exact in 
this limit. The cluster approximations give a much bet- 
ter agreement with Monte Carlo simulations, but it is 
not clear yet whether it becomes exact or not. However 
one can expect a tendency towards phase seperation due 
to 'attractive' interactions beween the cars. One stopped 
car can induce jams with a rather long lifetime since the 
restart probability of the first car is rather small. This 
tendency can already be observed in the cluster probabil- 
ities for p — 0.75 (see Fig. ^. The distributions become 
broader and especially at small densities the probabilities 
P(0, 1, 1) and P(l, 1, 1) which characterize the clustering 
are enhanced compared to the limit p —^ 0. 



5.4 Cluster Approximation vs. COMF 

In this section the results of the cluster approximation and 
COMF in special limits are compared. As already men- 
tioned, both methods yield the exact solution for v„iax — 
1. This is related to the fact that paradisical mean-field is 
exact in that case and both methods are able to identify 
all GoE states. For v^ax = 2 the situation is different. 
Here all three methods are only approximative. 

The deviations between COMF results and the Monte 
Carlo simulations can be understood by looking more closely 
at the microscopic structure of the stationary state. For 
p — there are three different stable structures: 00000, 
.1.1.1., and ..2..2....2...2.... The last structure comprises 
all configurations where all cars have velocity Vmax and 
at least Vmax empty cells in front. For fixed density there 
are many configurations which produce the same flux. In 
a large system these configurations are mostly made up 






10. Comparison of the 3-cluster probabilities 
ij+i,nj+2) obtained from computer simulations ( — ) 



P(0,0,1) 



Fig. 

and the 3-cluster approximation (■ ■ ■ ) for p 
occupation number of cell j. Note that P(1,0, 0) 
andP(l,l,0) = P(0, 1,1). 



of combinations of the three stable structure (with small 
transition regions). For p > 0, .1.1.1. is the most un- 
stable configuration. It tends to seperate into 00000 and 
2. .2. ...2. ..2 under fluctuations of the headways. COMF is 
not able to account for these fluctuations since headways 
are assumed to be independent. Therefore the weight of 
the conflgurations .1.1.1. - and therefore gi - is overesti- 
mated. In the 3-cluster approach, on the other hand, the 
headways are not independent. Therefore it is able to iden- 
tify the dominating local structure for p — > correctly. 

COMF is not able to describe the 'attractive' part of 
the interaction properly. T his can also be seen in the jam- 
size distribution (see Sec. 4^). COMF always predicts a 
strictly monotonous distribution ( |38[ ) , whereas the cluster 
approximation in principle is able to describe situations 
where the jam-size distribution has a maximum at some 
(small) value of n > 0. 



6 Summary 

Although cellular automata are designed for efficient com- 
puter simulation studies, an analytical description is possi- 
ble, although difficult. We have presented here four differ- 
ent methods which can be applied to CA models of traf- 
fic ffow. The first approach, a simple mean-field theory 
for cell occupation numbers, is insufficient since the im- 
portant correlations between neighbouring cells (e.g. the 
particle-hole attraction) are neglected. We therefore sug- 
gested three different improved mean-field theories. These 
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might be more reasonable to consider p = and p — I 
as fixed points in order to understand the behaviour for 
fixed Vmax ■ For p — there is a continous phase transition 
from laminar flow to a congested phase (see [^l||2^ and 
references therein) at c = l/{vmax + 1). This transition 
turns into a crossover at < p < 1 |£l|,^ . It would 
be desirable to investigate the behaviour of this transition 
close to p = 1 in more detail. 

The methods presented here can also be used for other 
CA models, e.g. variants of the NaSch model |22,E3,Eq,pj, 
pq,E9|. The investigation of the NaSch model has led to a 
better understanding of their advantages and limitations 
so that it is easier to choose the approach most suitable 
for a given problem. 
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Fig. 11. Same as Fig. M, but for p = 0.75. Here only the 
results of the computer simulations are shown. 



approaches take into account certain correlations between 
the cells. The simplest method is the so-called 'paradisical 
mean-field' theory which is based on the observation that 
certain configurations (Garden of Eden states) can never 
be generated by the dynamical rules due to the use of par- 
allel dynamics. The cluster approximation, on the other 
hand, treats clusters of a certain size exactly and cou- 
ples them in a self-consistent way. Therefore short-ranged 
correlations are taken into account properly. In contrast, 
car-oriented mean-field theory is a true mean-field theory, 
but here one uses a different dynamical variable, namely 
the distance between consecutive cars. In that way, certain 
correlations between cells are taken into account. 

All three improved MP theories become exact for Vmax = 
1. For larger values of Vmax they are just approximations. 
In principle, the cluster approximation and COMF (in 
combination with a cluster approach) can be improved 
systematically. This is, however, very cumbersome. 

An interesting observation is that the qualitity of the 
approximation depends strongly on the value of p. This in- 
dicates that the physics changes with p, contrary to com- 
mon believe. Evidence for this scenario comes from the 
microscopic structure of the stationary state. In the limit 
p — > it is dominated by repulsive interactions which tend 
to aligne the vehicles at a headway of at least Vmax cells. 
On the other hand, for p — > 1 there is a tendency towards 
phase seperation. Cars which had to stop due to a fluc- 
tuation will stand for a rather long time and thus lead to 
the creation of jams. In the deterministic case this even 
leads to the existence of metastable states. 

In 1^ it has been suggested that the behaviour of 
the NaSch model is governed by two fixed points, namely 
p ~ and Vmax ~ oo. Our investigations show that it 
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